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Parasitic two-level tunneling systems originating from structural material defects affect the func¬ 
tionality of various microfabricated devices by acting as a source of noise. In particular, super¬ 
conducting quantum bits may be sensitive to even single defects when these reside in the tunnel 
barrier of the qubit’s Josephson junctions, and this can be exploited to observe and manipulate the 
quantum states of individual tunneling systems. 

Here, we detect and fully characterize a system of two strongly interacting defects using a novel 
technique for high-resolution spectroscopy. Mutual defect coupling has been conjectured to explain 
various anomalies of glasses, and was recently suggested as the origin of low frequency noise in 
superconducting devices. Our study provides conclusive evidence of defect interactions with full 
access to the individual constituents, demonstrating the potential of superconducting qubits for 
studying material defects. All our observations are consistent with the assumption that defects are 
generated by atomic tunneling. 


INTRODUCTION 


Superconducting quantum bitd - 1 have recently achieved 
a breakthrough by demonstrating excellent gate fidelities 
and long coherence times in a fully scalable architecture,^ 
placing the realization of an integrated quantum comput¬ 
ing chip within reach. The solid-state approach however 
bears the burden that the material of the quantum device 
itself may host parasitic defects that give rise to two-level 
systems (TLS) acting as a sparse decohering bath. 

First signatures of coherent TLS in phase qubits were 
found in spectroscopy data, where observed avoided level 
cro ssing s manifest the defects’ two-level quantum charac¬ 
ter ! 3 10 Often, these defects show longer coherence times 
than the qubit itself 5 ^ and thus might be useful as quan¬ 
tum memorieJ^ and resources for quantum algorithms . 7 8 
Phase qubits were employed in several attempts to iden¬ 
tify the physical origin of those TLS, for example by ob¬ 
taining statistics on fre que ncies and coupling strengths , 2 
estimating their density,® measuring the temperature de¬ 
pendence of their coherence timesj^ or verifying theoreti¬ 
cal models describing their origin . 7 

The possibility of a direct interaction between TLS has 
been invoked in the past to explain the line width broad¬ 
ening and spectral diffusion of ultrasonically excited 
ensembles of TLS in glasse^ 1 ^ as well as var ious other 
low-temperature properties of disordered solids. 13 ^ TLS 
are furthermore a widely accepted model to explain noise 
in superconducting circuits, and mutual TLS coupling 
was recently suggested as the origin of the low frequency 
noise observed in microwave resonators.^ 


Here, we report the first clear experimental evidence 
of two coherently interacting two-level systems which are 




Figure 1. Using a superconducting qubit to access de¬ 
fects in Josephson junctions, a Illustration of the double¬ 
well potential for an atomic tunneling system. Tunneling en¬ 
ergy A and asymmetry energy 5 determine the level splitting 
E. b Sketch of the sample holder. To control the strain, the 
qubit chip is bent by applying a voltage V p to the stacked piezo 
actuator, c Schematic of the phase qubit including manipula¬ 
tion and measurement circuitry. The Josephson junction (JJ) 
tunnel barrier is sketched as a disordered insulator hosting 
TLS defects, here pictured as atoms tunneling between two 
metastable positions, with the arrows illustrating their elec¬ 
tric dipole moment. 


residing in the tunnel barrier of a Josephson junction. 
The data are obtained with a new technique for high- 
resolution defect spectroscopy that exploits the tunabil- 
ity of TLS by mechanical strain and their strong coupling 
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to a superconducting qubit. To characterize the coupled 
defect system in more detail, we extend upon this tech¬ 
nique and perform coherent two-photon spectroscopy that 
directly reveals the TLS’ coupling strengths and inde¬ 
pendent parameters. Interpretation of the measurement 
based on atomic tunneling systems fully accounts for all 
observations. 


RESULTS 

ATOMIC TUNNELING SYSTEMS 

To explain the microscopic origin of the TLS in su¬ 
perconducting electronics, several models® 20 have been 
proposed. However, all experimental results obtained so 
far, including the recent demonstration that the energy of 
the TLS is tunable by static mechanical strain, 1 ! are read¬ 
ily explained assuming that they originate from atomic 
tunneling systems. As in the well studied model describ¬ 
ing the low temperature thermal, dielectric, and acous¬ 
tic properties of disordered solids p- 2 *3 it is assumed that 
some atoms or small groups of atoms are able to tunnel 
between two energetically almost equivalent sites within 
the disordered oxide material of the device. These sys¬ 
tems give rise to two-level excitations in a wide energy 
range of up to the order E ~ /c#-lK or ~ h- 20 GHz. In 
bulk disordered solids, TLS are found in large numbers 
but, in contrast to their counterparts present in super¬ 
conducting qubits, cannot be addressed individually. 

According to the tunneling model, an atomic tunneling 
system is described as a particle in a double-well potential 
as shown in Fig. UK The energies of the two wells differ 
by the asymmetry 6 and the tunneling amplitude between 
them is denoted as A, resulting in a level s plitting of 
the two eigenstates given by E = \/ A 2 + s 2 .Tunneling 
systems couple to the environment predominantly by 
variation 5s of their asymmetry energy with 5s depending 
linearly on strain fields and, if the tunneling entity moves 
a non-zero charge, as well on electric field - the latter 
serving as an apparent explanation for the observed 
coupling of the TLS to the qubit circuit. A variation 5A 
of the tunneling amplitude induced by strain or external 
fields is generally believed to be negligible? 2 ^ We have 
recently verified the linear strain dependence of s and the 
corresponding hyperbolic variation of the energy splitting 
E by tracking individual TLS with a phase qubit while 
bending the chip circuit with a piezo actuator. 1 In our 
a setup, sketched in Fig. [TJd, an applied piezo voltage 
Up results in variable strain fields on the order of 10 _6 /U. 


DEFECT SPECTROSCOPY 

In this study, we detect and analyse TLS using a super¬ 
conducting qubit. These devices rely on Josephson junc¬ 
tions (JJ) as nonlinear circuit elements, which are realised 


as two superconducting films separated by a thin, insulat¬ 
ing tunneling barrier, consisting of a 2-3 nm thick struc¬ 
turally disordered layer of aluminum oxide. A sketch of 
the employed phase qubitP^ including measurement and 
manipulation circuitry is shown in Fig. [1J. The qubit’s 
level splitting and their population are controlled by ex¬ 
ternally applied flux bias and resonant microwave pulses, 
respectively, and a DC-SQUID is used for qubit readout. 

In order to trace the energies of individual TLS while 
applying strain to the qubit chip, we use a spectroscopy 
scheme based on the pulse protocol depicted in Fig. [2 ^l. 
The qubit is first biased at a frequency far away from the 
intended spectroscopy region and excited by a resonant 
microwave 7r-pulse. Applying appropriate flux bias, it is 
then tuned to the probing frequency fa where it resides 
for the holding time r. If at this frequency the qubit is in 
resonance with a certain TLS, the excitation is shared be¬ 
tween the systems.®! This results in coherent oscillations 
which effectively swap the quantum states of the two sys¬ 
tems at a frequency determined by their coupling strength 
as shown in Fig. It- A change 5P in the qubit excitation 
probability, measured after the interaction time r, thus 
reveals the presence of a TLS. We chose r to be about 
half the qubit’s life and coherence times T\ and T 2 (here, 
both ^100 ns) to reach a compromise between the loss of 
signal due to qubit relaxation and the sensitivity needed 
for detecting also weakly coupled TLS. We repeat this 
procedure for a range of probing frequencies fa and vary 
the mechanical strain applied to the sample. 

By plotting the change in qubit population as a func¬ 
tion of both applied strain and probing frequency, we ob¬ 
tain defect spectra like the one shown in Fig. §i- Dark 
traces indicate the resonance frequencies of individual 
TLS, which are tuned by strain as expected for atomic 
tunneling systems. Some TLS have a tunneling energy 
A that falls within the frequency range accessible by the 
qubit (about 6.5 - 10 GHz for our sample), while also 
their asymmetry energy 5 is tuned through zero in the 
investigated strain range. Accordingly, for those TLS we 
can clearly observe the hyperbolic strain-dependence of 
their resonance frequencies around minima given by A/h. 
We note that the distribution of TLS resonances changes 
completely once the sample is warmed to room temper¬ 
ature, see Supplementary Figure 1 for various examples. 
This can be explained by a modification of the atomic 
configuration changing the TLS environment locally, and 
also by an offset in the applied strain due to thermal 
dilatation of the sample fixture. As long as the tempera¬ 
ture is kept below about 10 - 20 K, the properties of the 
majority of TLS remain constant over several months of 
measurements. 

For TLS that are strongly coupled to the qubit, the 
chosen interaction time r may exceed the duration of 
one swap operation. Since the latter also depends on the 
detuning, a fringe-like interference pattern occurs around 
the traces, see Supplementary Figure 2. More details on 
this effect and other artefacts in such defect spectra are 
discussed in Supplementary Note 1. 
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Figure 2. Defect spectroscopy, a Pulse sequence used to detect TLS via resonant interaction with the qubit. The qubit is 
excited by a tt pulse, tuned to varying probe frequencies using different flux-pulse amplitudes A, and its population is measured 
after an interaction time r. b Qubit population probability, measured for two different probing frequencies in dependence of 
the interaction time r. The blue curve shows purely exponential energy relaxation for the isolated qubit, the red curve displays 
oscillations due to a strongly coupled and coherent TLS. We take the difference between the curves, measured at a fixed r as 
indicated, for the defect signal SP. c Defect signal 5P in dependence of the probing frequency for a fixed r. Individual TLS 
appear as pronounced dips, d Strain-dependence of TLS resonance frequencies, appearing as dark traces in 5P indicating a 
reduction in qubit population due to its resonant coupling to a TLS. Mutual TLS interactions are observed as random switching 
of the TLS resonance frequency (arrows), avoided level-crossings and non-hyperbolic traces (box). The cross-section at V p =-5 V 
(vertical line) is shown in c. 


MUTUALLY COUPLED TLS 

In the defect spectroscopy example shown in Fig. ir 
two different effects of interactions between TLS can be 
identified: i) The resonance frequency of certain TLS 
is observed to switch between two values in a random, 
telegraph-signal like fashion, with a switching frequency 
that depends on the applied strain. This effect is ex¬ 
plained by assuming that the observed TLS couples non- 
resonantly to an incoherent defect which fluctuates be¬ 
tween its positions, ii) Much less frequently observed are 
non-hyperbolic strain responses in the TLS resonance fre¬ 
quency as well as level splittings that are the typical sig¬ 
nature of two resonantly coupled coherent quantum sys¬ 
tems. 

In the remainder of this article, we discuss a partic¬ 
ularly clear manifestation of a system of two coupled 
TLS, whose defect spectroscopy signature revealed the 
5-shaped trace with avoided level crossings shown in 
Fig- i To explain this spectroscopic feature, we construct 
a model of two interacting defects denoted as ”TLS1” and 
”TLS2” (see also inset to Fig. |3|. Qualitatively speak¬ 
ing, we observe that TLS1 couples only weakly to the 


applied external strain so that the step-like increase of 
its resonance frequency is almost exclusively due to its 
non-resonant coupling to TLS2. In the region shown, the 
double-well potential of TLS2 is strain-tuned through its 
symmetry point ^(fp) = 0 at V p « — 14U such that the 
probability density of its ground state is shifted gradually 
from one potential well to the other. The accompanied 
shift in atomic positions is mediated via internal strain 
or electric field to TLS1, which responds by modifying its 
energy splitting, here from 6 to 6.8 GHz. Eventually, the 
applied strain tunes TLS2 into resonance with TLS1 and 
their coherent interaction gives rise to the associated and 
observed level repulsions. 

In the following, we outline our theory of the coupled 
defect system and show how it can be fully characterised 
by analysing strain-spectroscopic data. The Hamiltonian 
of a single TLS is written as 

Hi ~£i(y p )(j z ^ T — i =F/^(Fp) d Zj i, (1) 

where Ai is the tunneling energy and the asymmetry en¬ 
ergy £i(V p) depends linearly on external strain, i.e. volt¬ 
age Up of the piezo drive. Here and in the following, we 
use the tilde to distinguish operators such as the Pauli- 
matrices &j in the eigenbasis from those in the localized 
basis <jj. The energy splitting in the diagonal basis is 
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Figure 3. Spectroscopic signature of two mutually in¬ 
teracting TLS. The observed ^-shaped feature with avoided 
level crossings (encircled) is characteristic for a TLS (TLS1) 
that is coherently coupled to a second defect (TLS2). The 
inset shows their spectrum calculated from theory, closely re¬ 
producing the data obtained between 6 and 7 GHz (dashed 
box). Here, TLS2 is strain-tuned through its symmetry point 
£ 2=0 and changes the location of the ground state in its 
double-well potential (insets), while TLS1 is only weakly in¬ 
fluenced by strain. 

Ei(Vp) = y/si(V p ) + A?. The transformation to the di¬ 
agonal basis corresponds to a rotation about the angle 
£i(Vp), defined by tan^ = Ai/si with sin & = Ai/Ei and 
cos = Si/Ei. For example, the operators cr Zj i, whose 
eigenvalues identify the particle positions, transform as 

cr z? i -)► cos & d Zji + sin & cf Xji . (2) 

We write the Hamiltionian of two coupled TLS as H T = 
Hi + H .2 + i?i 2 , with the interaction term 

H \2 = - g cr z i cr Zj 2 - (3) 

Within the tunneling model, the defect’s mutual cou¬ 
pling parameter g comprises their electric dipole inter¬ 
action as well as a strain-mediated elastic contribution. 
Our spectroscopic data does not allow us to distinguish 
between these coupling mechanisms. However, since TLS 
interact with the qubit only electrically, we can determine 
the projections of the TLS’s electric dipole moments onto 
the field in the qubit junction independently. Details of 
this analysis are included in the Supplementary Note 2. 

In the diagonal basis, the interaction Eq. (J3| consists 
of four terms, which are combinations of <r z and cr x for 
each TLS (see Supplementary Note 3). However, for 
explaining the observed S-shaped signal (Fig. [3|, only 
two terms are relevant. The energy shift of TLS1, i.e. 
the amplitude of the ”5”, is given by the longitudinal 
coupling component g\\ = 2g cos cos £2 which stems 
from the term oc d z qd Zj 2 - The transversal component 



Figure 4. Exploring the energy level structure for 
the coupled defect system, a Level scheme of the reso¬ 
nantly coupled defect system for an applied strain near the 
right avoided level crossing of Fig. [ 3 ] at V p ~ 7V. Energy 
splitting and offset of the entangled states |=b) are determined 
by transversal and longitudinal coupling strengths g± and g \\, 
respectively, b Pulse sequence used to map out the complete 
energy level structure. For the upper panel in c, the qubit was 
measured after the first swap operation at varying frequencies 
in order to calibrate the entangled states energies. The lower 
panel was obtained using the complete sequence: preparation 
of one of the entangled states |—) or |+) in a first swap, fol¬ 
lowed by a second qubit excitation and variation of the second 
swap frequency to reveal the transitions between the entan¬ 
gled states to the fully excited state |ee). The colour maps 
were chosen according to the coloured arrows which indicate 
the corresponding transitions in a. Lines in c are calculated 
from theory. 


results in exchange coupling oc g ± 0-^,1 and defines 

the size of the level repulsion g_\_ = g sin £1 sin £ 2 - The 
remaining two parts of # 12 , oc d Zj id X; 2 and oc d Xj id z? 2 , 
yield only minor energy shifts and can be neglected 
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to first order. However, in our numerical fits to the 
spectrum we take the full interaction Eq. ©> into account. 


TWO-PHOTON SWAP SPECTROSCOPY 

To fully explore the coherently coupled defect system 
with even higher precision, we focus on the region near the 
right anti-crossing where both TLS are strain-tuned into 
resonance. This results in the four-level energy spectrum 
illustrated in Figure [4^i, which we map out by perform¬ 
ing microwave swap spectroscopy on the system prepared 
in different entangled states. For this, we follow the se¬ 
quence sketched in Fig. [4 |d. The qubit is first prepared 
in its excited state and then tuned to either the lower 
or upper branch of the avoided level crossing, realizing a 
swap operation with the entangled state |—) or |+), re¬ 
spectively (see §o- Directly afterwards, the JJ-qubit is 
again excited and tuned through a lower frequency range 
in order to find the transition that brings the system of 
coupled TLS to the fully excited state |ee). 

Data in the upper panel of Fig. © were obtained by 
measuring the qubit after the first swap operation, thus 
indicating the avoided level crossing as in our usual defect 
spectroscopy protocol in order to calibrate the prepara¬ 
tion of the chosen entangled state. The complete sequence 
results in the data shown in the lower panel of Fig. 
where we plot the difference in qubit population between 
two experiments in which the TLS system was prepared 
in either one of the two entangled state by adjusting the 
frequency of the first swap pulse. This experiment clearly 
reveals the transition to the fully excited state in excellent 
agreement with theory (dashed lines in Fig. [ 4 ]:). More¬ 
over, we directly obtain the energies E\ and E 2 of the two 
unperturbed TLS as well as the longitudinal and transver¬ 
sal coupling strengths as indicated in Fig. [4^. A detailed 
analysis of this experiment is contained in Supplementary 
Note 4. 

The spectrum calculated from our theoretical model, 
shown in the inset to Fig. [3j allows one to fit all sys¬ 
tem parameters and reproduces the data with high accu¬ 
racy (see Supplementary Fig. 3). The inter-TLS coupling 
strength g = —872 MHz is calculated from its compo¬ 
nents \g±\ = 155 MHz and g\\ = —428 MHz, obtained in 
swap spectroscopy for the TLS tuned into resonance at 
V p = 7 V as explained above. 

Together with the energy splitting Ei of the individual 
TLS at their resonance, we can determine their strain- 
dependent mixing angles ^ to fully characterize the sys¬ 
tem. The obtained TLS parameters including their de¬ 
pendence on applied piezo voltage are summarized in Ta¬ 
ble [TJ more details on this evaluation are given in Supple¬ 
mentary Note 3. 

A final remark concerns the probability of finding 
two TLS spaced closely enough to expect an interaction 
of similar strength as in our experiment. Considering 
only electric dipole interaction and assuming for both 



A 

£i(V p) @£2=0 

TLS1 

5.47 GHz 

3.18 GHz - 4 MHz/V 

TLS2 

1.3 GHz 

295 MHz/V 


Table I. Measured parameters of the two coupled TLS. 

The asymmetry energies Ei are given at the symmetry point 
of TLS2 (V p = -14.05 V). 

TLS an electric dipole moment of d ~ 1 eA, which is 
consistent with the results of our work and agrees with 
recent observations and theory EESHI 20 we can make 
a rough estimate of the maximal distance between the 
two interacting TLS to be on the order of 5 nm (see 
Supplementary Note 4). When distributing the ^ 50 
TLS visible in the accessible frequency range of about 
6 to 9 GHz evenly onto the area of the 1 fim 2 large 
qubit junction, on average each TLS occupies an area 
of about 140x140 nm 2 . Thus, although one may expect 
an increased TLS density at interfaces, observing two 
coherently interacting TLS in Josephson junctions is 
indeed a rare case. 


DISCUSSION 

The experimental techniques presented here provide a 
novel spectroscopic view onto the bath of sparse mate¬ 
rial defects by accessing the quantum states of individual 
TLS and small coupled systems with a superconducting 
qubit while strain-tuning their internal degree of free¬ 
dom. This lays the ground for further experiments to 
clarify the microscopic origin of TLS, which is vital for 
the advancement of various kinds of nanofabricated de¬ 
vices whose functionality is hampered by defects. So far, 
the atomic tunneling model readily explains all effects 
observed with TLS in tunnel junctions and, in partic¬ 
ular, the here demonstrated strong coupling of TLS to 
mechanical strain and their interaction with both coher¬ 
ent as well as randomly fluctuating defects. Our results 
open way to detailed testing of the 50-year old tunneling 
model on the basis of individual TLS, for example by per¬ 
forming defect spectroscopy for statistical analyses of the 
TLS distribution and by studying the strain dependence 
of TLS coherence. 


I. METHODS: SUPERCONDUCTING QUBIT 
SAMPLE 

The phase qubit sample 25 used in this work was fab¬ 
ricated in the group of J. M. Martinis at University of 
California, Santa Barbara (UCSB). The qubit junction 
had an area of about 1/im 2 , fabricated using aluminum 
as electrode material and its thermally grown oxide as a 
tunnel barrier. All data have been obtained at a sam¬ 
ple temperature of about 35 mK. The mechanical strain 









6 


was controlled by bending the sample chip with a piezo 
transducer. 1 
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SUPPLEMENTARY NOTE 1: DEFECT 
SPECTROSCOPY 

Our defect spectroscopy protocol relies on detecting 
the transition of a single microwave photon, initially 
stored in the qubit, via resonant interaction to individual 
TLS. In comparison to our previous method, 1 where a 
long microwave pulse was saturating the qubit transition, 
here we use a 7r-pulse to maximize the qubit population 
probability and thus increases the signal-to-noise ratio. 
In principle, this method equally detects strongly coupled 
circuit or environmental resonances, which can however 
be distinguished by their independence on strain. We 
did not observe any strain dependence of the qubit 
parameters nor presumed circuit resonances. We note 
that a similar, independently developed technique of TLS 
swap spectroscopy was used in RefP to obtain statistics 
on TLS coherence and coupling strength, but without 
mechanical strain control. 

Observation of coupled TLS systems 

We performed defect spectroscopy on the same sample 
and setup in 12 different cool-downs over a course of two 
years. While the qubit parameters remained constant, 
the distribution of TLS resonances changed completely. 
This is not surprising given the observed strong coupling 
of TLS to the local mechanical strain, which is expected 
to vary between cool-downs due to small reconfigurations 
of atomic positions. However, changes in defect spec¬ 
troscopy can also be explained by offsets in the exter¬ 
nally applied strain due to thermal dilatation affecting 
the sample holder. 

5-shaped signatures of coherently coupled defects, 
shown in Supplementary Figure [l] were observed in two 
different cool-downs. 

Supplementary Figure & was obtained in the exper¬ 
imental run in which we characterized the coupled de¬ 
fect system as discussed in the main text. About one 
year later, the same sample revealed a similar system of 
two coupled defects shown in the red box in Efc having 
comparable longitudinal coupling strength but different 
resonance frequencies and opposite symmetry. It may 
well be that both observations involve the same tunnel¬ 
ing systems which however experience slightly different 
local potential configurations. 

In Supplementary Figure [ljo, one can notice that the 
frequency of the coupled defect system switched while the 
spectrum was measured around a voltage of about -25 V 
(the dashed red arrow indicates the shift in frequency). 
Such spontaneous changes in the resonance frequencies 
of some TLS are observed frequently and occur most 
often shortly after the sample was cooled to milikelvin 
temperatures. This is readily explained by a coupling of 
the observed coherent TLS to incoherent TLS, also called 
” two-level fluctuators”, which may be of same physical 
origin but have a small rate of tunneling between their 
potential wells. If such a fluctuator is strongly asymmet¬ 
ric, it may be found in its higher energy state soon after 


the cool-down while it would remain in its ground state 
once tunneling occurred. Experimentally, we found that 
the number of spontaneously switching TLS is reduced 
by repeatedly cycling the applied strain through the 
complete range, hereby ” annealing” the sample in the 
sense that fluctuators are stimulated to tunnel to their 
more stable ground state. The data shown in Fig. 2 d of 
the main text has been measured in the same cool-down 
as [l] b, but 3 weeks later. One can see that some TLS 
changed their properties in between measurements, while 
others remained stable. More examples of mutual defect 
coupling are seen in Supplementary Figure showing 
telegraphic switching of a TLS’ resonance frequency 
(arrow) and an avoided level crossing (circle). 

Effects and experimental artefacts 

In the following, we discuss some artefacts that may be 
observed with this spectroscopy technique. 

Fringes or shadow-like replica of some TLS traces, such 
as shown in Supplementary Figure UK may appear for 
TLS whose coupling strength to the qubit exceeds h/2r. 
In this case, the chosen interaction time r exceeds the 
duration of one swap operation such that the excitation 
is transferred back to the qubit, and this lowers the de¬ 
fect signal SP (see also Fig. 2b in main article). Since the 
frequency of swap oscillations increases with the detun¬ 
ing between qubit and TLS, a minimal qubit population 
(maximal signal |5P|) thus occurs only in the vicinity 
of the exact resonance as shown in Supplementary Fig¬ 
ures 03 and c. 

Fringes may also occur due to residual entanglement 
between the qubit and TLS. An example is visible in 
Supplementary Figure [lfc in a range of V p « 40...50 V 
at low frequencies. In this measurement, the qubit was 
biased at a frequency of 6.4 GHz, so that already during 
excitation it was near resonance with a strongly coupled 
TLS. Accordingly, the 7r-pulse prepares some entangled 
state between the systems, whose phase is then modified 
in dependence of the subsequent qubit detuning during 
the interaction time r. 

Qubit drift due to uncontrolled slow changes in its bias 
flux leads to deviations in the defect signal’s frequency de¬ 
pendence. For Supplementary Figure [T]^, in total about 
140,000 data points were measured, of which each was 
averaged 1000 times at a repetition rate of about 500 
Hz, resulting in a total measurement duration of ~ 3.5 
days. An example of qubit drift can be seen in Supple¬ 
mentary Figure [TJ^ around a piezo voltage of about 7 V, 
which was measured while feeding liquid Helium to the 
dilution refrigerator. This resulted in a small change of 
the qubit flux bias because of temperature variations in 
its bias line filters that are installed at the lK-pot of 
the cryostat. Qubit drift can be distinguished from TLS 
resonance frequency fluctuations since all TLS signals are 
shifted equally. We note that the measurement time could 
be easily reduced by at least one order of magnitude if 
one employs a dispersive qubit readout methocP® instead 
of the DC-SQUID switching-current measurements em- 


ployed hereP 

SUPPLEMENTARY NOTE 2: THE QUBIT-TLS 
COUPLING 

Assuming the TLSs interact with their environment via 
their electric dipole moment with the qubit circuit, the 
corresponding coupling operator is o z , whose expectation 
values identify the position of the particle. The interac¬ 
tion of an individual TLS with the electric fields inside 
the qubit’s junction is then described as 

H iq = ^Vif x a Zj i , (4) 

where Vi is the electric dipole interaction strength and 
we used the fact that for a phase qubit in its eigenbasis, 
the electric field operator is ocr*, P Rewriting Supplemen¬ 
tary Eq. @ in the TLS eigenbasis, we retain two terms, 
of which the transversal coupling term Vi : _\_f x a x ^ with 
Vi^± = Vi sin describes the resonant interaction between 
the JJ qubit and the TLS i. It is this term in particular 
which causes the exchange of excitations between the two 
systems. 

We obtained the electric coupling strength between 
qubit and TLS1, vi t ±, by fitting time-domain oscillations 
of the initially excited qubit tuned into resonance with 
TLS1, far away from the TLS1-TLS2 anti-crossing, where 
uq_L is independent of the state of TLS2. We performed 
the same experiment also with both, lower and upper 
branches, in the vicinity of the TLS1-TLS2 anti-crossing 
(the pulse sequence is shown on the inset in Supplemen¬ 
tary Figure | 4 ]). The result is plotted in Supplementary 
Figure green for the lower and red for the upper 
branch. The panels in Supplementary Figure and c 
show the observed oscillations in the qubit population. 

Since TLS2 is very weakly coupled to the qubit, the os¬ 
cillation frequency decreases quickly with voltage when 
following a branch turning into TLS2. However, it is 
worth noting that the crossing point of the two curves 
is not at the voltage of 7 V where TLS1 and TLS2 
are in exact resonance. This would be expected if the 
qubit-TLS2 coupling was zero, as indicated by the dashed 
curves. A non-zero v 2 ,± explains this shift, because the 
coupling between the qubit and the |=b) states scales as 
{ v i,± 4= U2,±)/V / 2. Then, the theory shows a very good 
agreement with the experiment and we obtain v^± = 15.4 
MHz and v 2 ^± = 3.0 MHz at an applied strain where the 
TLS are in resonance E\ = E 2 . 

With the already known £1 and £ 2 , Vi = 18.7 MHz 
and V 2 = 14 MHz can be calculated from Supplemen¬ 
tary Equation 0. It is worth mentioning that the very 
small v 2 ,± is not attributed to a smaller dipole moment 
in direction of the qubit electric field, but can be simply 
explained by a strong asymmetry A 2 £i resulting in a 
very small dipole moment of the TLS eigenstates. This 
analysis is restricted to the assumption that both TLSs 
interact purely electrically. 


SUPPLEMENTARY NOTE 3 : ANALYSIS OF 
THE FULL HAMILTONIAN 

We write the full Hamiltonian describing the JJ-qubit 
coupled to two interacting TLSs as 

2 2 

H = H q + Y,ni + Y, H * + H ' 2 ’ ( 5 ) 

i=1 i=l 

with the qubit Hamiltonian H q = | E q f z , where E q is 
the level-splitting of the two lowest qubit states and f z 
is a Pauli-matrix describing the qubit eigenstates. Hi is 
the Hamiltonian of the i th TLS, while H iq describes its 
interaction with the qubit, and H\ 2 the TLSs’ mutual 
interaction as described in the main manuscript. 

Diagonalizing the Hamiltonian in Supplementary 
Eq. ([ 5 ]) introduces several new terms. The resulting 
Hamiltonian, if, can be significantly simplified by ignor¬ 
ing all coupling terms of the form oc a x a z and oc cr z cr x . 
They represent a coupling where one partner changes its 
state depending on the instantaneous state of the other 
subsystem and contribute only as small energy offsets. 
Neglecting these minor energy shifts, we can write the 
significant parts of the full Hamiltonian as 

H = Eqf z + E\g z ,\ +E 2 cr z , 2 + ^||^,i^,2 

E 9-L&x,i&x,2 + E ^ 2 ,a_^x,2^~x )5 ( b ) 

with 

Ei = s[a + 4 2 > 

= Vising , v 2 ,x = ^2 sin £2 , 

5|| = 5 cos £1 cos £2 , g±. =9 sin £1 sin & , 

cos €i=ei/Ei , sm£ i = A i /E i (7) 

Here the relationship between the longitudinal and 
transversal coupling factors g\\ and g_ i_, which are eas¬ 
ily identifiable in experiment, and the TLS mixing an¬ 
gles & becomes clear. In our case, the ^ dependence 
of g± cannot be observed because when TLS2 is de¬ 
tuned from TLS1, the perpendicular coupling oc o Xi \o x ^ 2 
becomes irrelevant We can interpret the term cos ^ as 
representing the expectation value of the position op¬ 
erator in the double-well potential (cr^^), so that g\\ = 
9 (&z,i) (&z,2)- As explained in the main manuscript, 
when tuning through the TLS2 symmetry point, (cr Z:2 ) 
changes its sign, corresponding to reversing the direction 
of TLS2’s dipole moment. Due to the strong dipole-dipole 
interaction, TLS1 adjusts its energy splitting accordingly, 
and this results in the observed 5-shaped signature. 

Additionally, the angles & also change the qubit-TLS 
coupling strength Vi r ± which determines the visibility of 
TLS in our defect spectroscopy. The coupling strength 
Vi ,_l is maximal when the TLS is at its symmetry point 
Si = 0 , and it is strongly suppressed for largely asymmet¬ 
ric TLSs with Ai Si. These considerations allow us to 
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obtain a close fit of the theoretical model to the experi¬ 
mental data, with the results displayed in Supplementary 
Figure [ 3 ] and table I, respectively. 

During our measurements, we occasionally observe 
abrupt jumps of individual TLS’ resonance frequencies on 
a time-scale of days or weeks. From defect spectroscopic 
data obtained before and after such events, we found that 
the defects experience changes of both their asymmetry 
and tunneling energies, while the strain dependence of 5 
remains very similar. A sudden change in TLS param¬ 
eters also occurred between measurements of the defect 
spectrum and that of the fully excited TLS1-TLS2 system 
(see Appendix B), by which the TLS1 energy shifted by 
approximately 100 MHz towards lower frequencies. Since 
after the jump, the theoretical spectrum is found to still 
coincide with the anti-crossing on the left side of the “S” 
(data not shown), we can conclude that A 2 did not change 
appreciable in comparison to £2 at the TLS1-TLS2 anti¬ 
crossing. 


SUPPLEMENTARY NOTE 4: MEASURING THE 
ENERGY OF THE FULLY EXCITED TLS1-TLS2 
SYSTEM 


The aim of the experiment discussed here is to measure 
the transition energies of the coupled TLS1-TLS2 system 
around the anti-crossing at F p = 7 V on the right side of 
the “S”-shaped signature shown in Supplementary Fig¬ 
ure [ 3 ] The energies of the four levels \gg), |=b) and |ee) 
are illustrated in Fig. 4 a. While the splitting size equals 
the energetic shift of the anti-crossing from E ee / 2, 
which corresponds to the unperturbed TLS energies E\ 
and E 2l yields gy A^Fig. 4 b visualizes the pulse sequence 
for this experiment. The upper image of Fig. 4 c shows 
a zoom of the right anti-crossing with the two branches. 
After exciting the qubit with a 7r-pulse, the excitation 
is swapped to one of the branches by tuning the qubit 
for the swap time to the corresponding energies around 
E+ or E-. Afterwards, the qubit is excited a second 
time and tuned to lower frequencies in order to search 
for the transition that fully excites the TLS system. This 
occurs at energies around E ee — E+ or E ee — E _, respec¬ 
tively, and the qubit will loose its excitation yielding an 
additional dark trace in the spectrum image (data not 
shown). The two plots obtained in this way, one for each 
branch into which the first excitation was swapped, are 
subtracted from each other yielding peaks or dips at the 
relevant energies and zero for the background (Fig. 4 d, 
lower plot). For better visibility, a color-map has been 
chosen such that peaks in Fig. 4 c, which correspond to 
first exciting the lower branch, appear in green and the 
dips, arising if the qubit excitation was first swapped to 
the upper branch, show up red. The coupling constants 
between the two TLSs, and also the TLS energies in the 
uncoupled case, can now be directly extracted from the 


data: 


2E\ — 2E 2 — E+ + (E ee — E+) — E + (E ee — A_) 
g± = A + — A = (E ee — E_) — (E ee — Efa) 

2< 711 = E+— ( E ee - A_) = E_ - (A ee - E+). ( 8 ) 

With the known ratio of the tunneling and asymmetry 
energies of TLS2, £ 2 ? g and £1 (Supplementary Eq. 0 ) 
can be calculated (see Tab. I). 


SUPPLEMENTARY NOTE 5: DISTANCE 
BETWEEN THE TLS 

In order to obtain a rough estimate of the distance 
between the coupled TLS, several assumptions have to be 
made. First of all, for the calculation we assume that the 
mutual TLS interaction is solely due to electrical dipole 
coupling, i.e. we neglect any contribution due to elastic 
coupling. The interaction strength can then be written 

aJ3 

2 = 4neo£r^ — 2 Mi||M 2 l|) > (9) 

where r is the relative position vector from TLS1 to TLS2 
and their electrical dipole moments fa = fa^ + fa f ± are 
decomposed into components parallel and perpendicular 
to r, respectively. 

We estimate the magnitude of the TLS’ electrical dipole 
moments by their components parallel to the electric field 
in the Josephson junction, which can be determined as 
described in Supplementary Note 3. The qubit-TLS cou¬ 
pling strength v^\\ is the product between the TLS dipole 
moment and the electrical field in the junction, 10 


0 Mi, \\ a [Ew 

2 ~\j2C 


( 10 ) 


where x ~ 2 nm is the thickness of the Junction’s tun¬ 
nel barrier, C = 850 fF the qubit’s capacitance and 
A 01 ~ h • 6.5 GHz the qubit’s energy level splitting at 
resonance with the TLS. Using the values of Table I, we 
obtain /i 1? y « 0.46 eA and /i 2 ,\\ ~ 0.38 eA. Finally, we 
assume yiq = (i^ y and use Supplementary Eq. ^ with 
e r = 10 (saphhire) to estimate the maximal distance be¬ 
tween the TLS, which results in r = 4.1 nm and r = 5.2 
nm, for horizontal and vertical arrangements of the TLSs, 
respectively. We note that this maximal vertical distance 
exceeds the thickness of the tunnel barrier. 
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strain / piezo voltage t/ p (V) 


Supplementary Figure 1 . Defect spectroscopy examples obtained on the same sample. Red boxes highlight coupled TLS 
systems, circles indicate avoided level crossings, a Data taken in April 2012, showing partly the signature of the coupled TLS 
system that was characterized in this work, b Data acquired soon after cooling down the sample in March 2013, about 3 weeks 
before taking the data in Fig. 2d of the main manuscript. A few TLS changed their properties, while most remained stable. 
The dashed arrow indicates a spontaneous change in TLS parameters, c Another example of an observed avoided level crossing 
(circle), and telegraphic switching (arrow). 
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interaction time z 0 -0.5 

5P, defect signal 


Supplementary Figure 2. Origin of fringes in defect spectra.a Defect spectroscopy example, showing the hyperbolic 
resonance frequency traces of three TLS. Fringes appear around the exact TLS resonance for strongly coupled TLS (middle 
and rightmost hyperbola), b and c Calculated population probability P(|l)) of an initially excited qubit, showing coherent 
oscillations due to coupling to a TLS in dependence of the interaction time r and detuning (left panels). The right panels 
show the defect signal 5P for a fixed interaction time r (cross-section along vertical white line in the left panel). For c, the 
TLS-qubit coupling strength is twice larger than for b, resulting in faster oscillation. For the more strongly coupled TLS, 
but same interaction time r, the signal 5P reaches minima at non-zero detuning, giving rise to the pronounced fringes or 
shadow-lines in defect spectroscopy. 
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Supplementary Figure 3. Comparison between experiment and theory, a Strain-spectroscopy data showing the 
characteristic ^-shaped signature of two coupled defects. The calculated spectrum is superimposed on the data, showing the 
excellent agreement with theory, b A plot of the whole theoretical spectrum using the values from Tab. I. The dashed hyperbola 
visualizes the TLS2 trace invisible in spectroscopy, while the dotted curve shows the energies of the fully excited system. 
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Supplementary Figure 4. Measuring the coupling strength between the qubit and each TLS. a Measured 
frequencies of oscillations in the qubit population for the excited qubit tuned to lower and upper branches of the avoided level 
crossing, respectively, which equals the coupling strength to the corresponding entangled states |+) and |—). The inset shows 
the used pulse sequence. It can be clearly seen that the measured coupling strengths (circles) deviate from the dashed curves 
which are calculated for the case of V 2 ,± = 0. However, one can account for the shift by introducing a non-zero qubit-TLS2 
coupling (solid lines), b and c: Observed coherent swap oscillations in P(|l)) (colour-coded) vs. strain and the interaction 
time. Their decay is dominated by qubit decoherence. 
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